Dynamic Magnetoelectric Effect in Ferromagnet|Superconductor Tunnel Junctions 
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We study the magnetization dynamics in a ferromagnet|insulator| superconductor tunnel junction and the as- 
sociated buildup of the electrical polarization. We show that for an open circuit, the induced voltage varies 
strongly and nonmonotonically with the precessional frequency, and can be enhanced significantly by the su- 
perconducting correlations. For frequencies much smaller or much larger than the superconducting gap, the 
voltage drops to zero, while when these two energy scales are comparable, the voltage is peaked at a value de- 
termined by the driving frequency. We comment on the potential utilization of the effect for the low-temperature 
spatially-resolved spectroscopy of magnetic dynamics. 
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The field of spintronics has evolved tremendously over the 
last decades, leading to important conceptual and technologi- 
cal advances in spin-based memories, sensing, and logic Q. 
The manipulation and detection of the spin degrees of free- 
dom, such as collective magnetization, lies at the heart of spin- 
tronics, with magnetic field, static or time-dependent, provid- 
ing a direct way to access it. It turns out, however, that the 
electric rather than magnetic control can be often preferred 
for spintronic manipulations |2], as the former can exert larger 
torques, act faster, and can be applied or detected with a finer 
spatiotemporal resolution (3j|4]|. 

There are two main routes for the electrical control of mag- 
netization dynamics: One relies on spin-polarized electrical 
currents, which couple to the magnetization via the spin- 
transfer torque 1 3 ], while the other, which is also a more recent 
development, is based on controlling magnetic anisotropics by 
applying voltage pulses [4 ] . Focussing on the former, the most 
basic implementation for creating spin-polarized currents is to 
pass unpolarized electrical current through a fixed ferromag- 
net. An alternative way for generating spin currents relies on 
the spin Hall effect, which leverages spin-orbit interaction in 
the material and does not require any ferromagnetic polariz- 
ers. In fact, this has been established as a primary tool both 
for manipulating and detecting the magnetization dynamics, 
transforming spin signals (spin currents) directly into electri- 
cal signals (Hall voltages) or vice versa |5|. The spin-Hall- 
induced voltage scales with the lateral dimension of the sam- 
ple, making itself extremely useful for larger devices [6]. The 
main drawback, however, is that it loses its utility when it 
comes to detecting local magnetization dynamics. 

In this Letter, we study voltage induced by magnetization 
dynamics in a circuit involving a driven magnet coupled to an 
s-wave superconductor via a weak tunnel barrier. It was previ- 
ously shown that magnetization dynamics can induce voltage 
in tunnel junctions with normal metals or static reference fer- 
romagnets EH), by the process of an adiabatic charge pump- 
ing. The resulting voltage is, however, small, compared to 
the driving frequency, at typical microwave powers. We show 
here that singularity associated with the quasiparticle density 
of states can significantly enhance such dynamically induced 
voltages at driving frequencies corresponding to the supercon- 



FIG. 1 . A schematic of the system. The ferromagnet is driven to 
precess with frequency co at angle 6 with respect to the z axis. The 
effective splitting in the rotating frame of reference is co cos 6 along 
the magnetization direction m within the ferromagnet. The super- 
conductor is subjected to the effective spin splitting co along the z 
axis, which shifts the quasiparticle bands in the rotating frame. A s is 
(twice) the s- wave superconducting gap and \i is the chemical poten- 
tial in equilibrium in the laboratory frame. The magnetic precession 
pumps charge current / through the tunnel barrier in the closed cir- 
cuit, or a voltage V is measured by a voltmeter in the open circuit. 



ducting gap. The underlying pumping behavior is thus neces- 
sarily nonadiabatic. 

The Bogoliubov-de Gennes Hamiltonian for our hybrid fer- 
romagnet/superconductor junction (see Fig.[T]) reads [9] 
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in a certain basis, where m(f) is the magnetization direction 
in the ferromagnet, A F (r) = A F Q(-x) is the magnetic ex- 
change field, A s (r) = A s ®(x) is the (real- valued) supercon- 
ducting pair potential, both written in terms of the Heaviside 
step function @(x), a = (<x x , cr y , cr z ) and r = (r x , r y , r z ) are 
Pauli matrices operating in spin and particle-hole (Nambu) 
subspaces, respectively, and V(r) is the total effective scalar 
potential acting on electrons. We suppose the ferromagnetic 
and superconducting regions are separated by an insulating 
barrier, such that V(r) is large near x « 0, where x stands for 
direction normal to the junction placed at x = 0. 



For a circular precession, we parametrize the magnetization 
direction as m(t) = [sin cos (cot), sin 6 sin (cot), cos 6], with 
6l> being the precession frequency and 6 the precession angle. 
The exchange splitting A^ is assumed to be much larger than 
both the superconducting gap A^ and the precession frequency 
co: A F » As , co (setting ft = 1 throughout), which is typically 
the case. For a steady precession, it is convenient to switch to 
a rotating frame of reference, where the ferromagnet is static. 
This is achieved by a time-dependent unitary transformation 
of the Hamiltonian: H^ s = U\t)H ms (t)U(t) - iU\t)d t U(t) = 
#Fi S (0) - cocr z /2, with U(t) = exp(-icr z cot/2) t8l[TUl In the 
ferromagnetic bulk, the Hamiltonian in the rotating frame be- 
comes H' F = H F (0)-(co/2)(T\\ cos 0, where cr\\ = m(0)cr is the 
spin projection on the magnetization direction ra(0). We have 
disregarded the component perpendicular to m(0), which is 
effectively suppressed for co/A F «c 1. On the superconduct- 
ing side, on the other hand, the spin splitting in the rotating 
frame is simply given by co, with the corresponding Hamil- 
tonian H' sc = //sc(0) - cocr z /2. The spectrum of the super- 
conductor in the rotating frame is thus shifted by the effective 
magnetic (Larmor) field, E { [\e) = ± ^e 2 + (A s /2) 2 - cocr/2, 
corresponding to the quasiparticle density of states 

Dy(E) = D —^±- " Cr/21 (2) 

^(E + cocr/2) 2 -(A s /2) 2 

Here, Z) is the Fermi-level density of states per spin projec- 
tion in the normal state, and <x = ± for spins up/down along 
the z axis. D^(E) diverges at energies E — > (±A S - cocr)/2. 
Since the spin-up and -down quasiparticle sub-bands in the 
superconductor are shifted by co, the gap closes in the ro- 
tating frame for co > A s . The spin-dependent spectrum on 
the ferromagnetic side is shifted too, but by co cos 6 instead 
of co. In the tunneling regime, we can assume that the fer- 
romagnetic and superconducting bulks are in their separate 
equilibria in the laboratory frame. This relies on the fact 
that the magnetization-dynamics induced pumping is expo- 
nentially weak, such that, in particular, it has essentially no 
effect on the self-consistent pairing potential A^ . 

We now rewrite Hamiltonian ([T]) in the tunneling approxi- 
mation, which is given in the second-quantized form by 



Ht = 2_j tk^crCp-kcrCS^cr + H.C. , 
k,q,cr 



(3) 



in the rotating frame. Here, tk, q ,o- = t is the tunneling ma- 
trix element, which is taken to be constant, for simplicity, and 
CF{S);k,cr are the electron annihilation operators in the ferro- 
magnet (superconductor), k and q label orbital quantum num- 
bers and cr spin projection on the magnetic direction. We do 
not expect a general spin- and k, qr-dependent tunneling to af- 
fect qualitative features of our final results, apart from mod- 
ifying parameters associated with the spin-dependent density 
of states in the ferromagnet. Hamiltonian ^ can be used to 
calculate the charge and spin currents flowing from one metal 
to the other across the tunnel barrier. 

We compute the out-of-equilibrium charge current using 
Fermi's golden rule for the transition probabilities involving 



all electron- and hole-like branches in the superconductor. In 
the tunneling regime, the superconductor can be effectively 
viewed as a simple semiconductor (TT] [121 (spin split in the 
rotating frame), with a singular quasiparticle density of states 
at its band edges, according to Eq. ^. The total current in- 
duced by the magnetic driving in the presence of a voltage V 
reads: 



/ =2ne\t\ 2 J] fdED^Df'\E) |<cr|cr'>| 2 
x[ff\E-eV)-ff\E)\, 



(4) 



where V is the electrochemical potential (applied and/or in- 
duced) of the ferromagnet relative to the superconducting 
condensate, e < is the electron's charge, Dp is the 
spin-<x Fermi-level density of states in the ferromagnet, and 
f F a \E) = {exp[(E + cocrcos(6)/2)/k B T] + l}" 1 , ff(E) = 
{exp [(E + coo-/2)/k B T] + l}" 1 are, respectively, the Fermi- 
Dirac distributions in the ferromagnet and superconductor 
at temperature T. For the spin matrix elements, we have 
|( C r|cr , )| 2 = cos 2 (0/2)cW +sin 2 (0/2>W, where & = -cr. We 
are now equipped to calculate the resultant current in the pres- 
ence of the magnetization dynamics, at an arbitrary tempera- 
ture. In contrast to band semiconductors, the gap A^ — > A^ (T) 
itself depends on T, closing at T c « 1.76A^(0) within the s- 
wave BCS model [12], while the Fermi level is pinned, in the 
laboratory frame, midgap by the superconducting condensate. 
We start by computing the charge current as a function of co 
and V at T = 0. This regime allows for an analytical evalua- 
tion of the current, as well as captures qualitative features that 
extrapolate to finite T. Writing / = £o-o-' I a(T >, we arrive at: 



I<ro> =I ^\ l + w' cos 0) ^[co(cr - v' cos 6)12 - V] 2 - 1 

X | Yj ®M<r - cr' cos6>)/2 - V + s] - 1 1 , (5) 
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where f ' - neAsi^tyDp'D^. Here, we are measuring all 
the energies in units of the zero-temperature gap, A^ (0)/2. In 
the following, we focus on the open-circuit case (see Fig. 1), 
so that the current flowing through the heterostructure is zero 
(/ = 0), while the voltage V induced by the magnetization 
dynamics is measured by a voltmeter. We are specifically in- 
terested in regimes where this voltage can be significantly en- 
hanced by the superconductor. In the tunneling regime, the 
Andreev processes are strongly suppressed, so that transport 
is governed by the excited quasiparticles. This means that for 
co < 1 (in units of A^ /2) and T = 0, the quasiparticles are 
gapped out and V = 0, as long as the microwave power is 
low enough so that multimagnon processes do not contribute. 
(Although absent from our model with circular precession, we 
expect the multimagnon processes to become relevant for gen- 
eral elliptical precession at sufficiently high microwave pow- 
ers.) 

In Fig. [2j we plot the open-circuit voltage V at T = 0, as 
a function of frequency co and precession angle 6. We find 
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FIG. 2. Dependence of the induced voltage V on frequency co and 
cos 6 at zero temperature. The thick black curve shows the activation 
frequency cdq - 2/(1 + cos#), while the red and blue curves corre- 
spond to the voltages V c \ [Eq. {7}] and V c2 [Eq. |9])] that signal the 
activation of new transport channels in the current /. In these plots, 
we set P = 2/3 and expressed V and co in units of A s /2. 



three different regimes in the dependence of V on co, corre- 
sponding to the activation of new spin channels in the tunnel- 
ing current §5§ with increasing co. The voltage depends on co 
nonmonotonically, reaching a maximum V c \ which, as shown 
later, can be as large as 1 . More specifically, we find that be- 
low a minimum frequency coq = 2/(1 + cos#) the voltage is 
always zero (assuming only positive frequencies, i.e., co > 0), 
while for co > coq there is a finite voltage drop in the system. 
The total current in this case is given initially by 1^ + 1^ for 
< 6 < n/2 and 7 TT + 7 U for 7r/2 < 6 < n. At certain higher 
frequencies (whose specific values are discussed below), the 
terms 1^ and In get activated for < 6 < n/2 and % and 
fy for n/2 < 6 < 7i. Successive activation of these differ- 
ent tunneling channels define voltage landscapes as shown in 
Fig.0 

Note that the symmetry of our system dictates that V(-co) = 
-V(co) as well as V(n - 6) = -V(6), allowing us to henceforth 
restrict our discussion to 6 e (0, n/2] and co > 0. While it is 
possible to extract analytical expressions for the voltage (and 
the activation frequencies) as a function of frequency at T = 0, 
these are too long and unilluminating. Instead, we will derive 
approximate expressions for the induced voltage in different 
frequency ranges. 

At frequencies co above coq but still sufficiently low such 
that \co cos 2 (6/2) ± V\ > 1 and |^sin 2 (0/2) ± V\ < 1 [where 
V = V(co, 6, P) needs to be solved self-consistently for, in the 
open circuit], the voltage V increases monotonically with co, 
as shown in Fig.[2j given approximatively [as an expansion in 



(1 - P)] by the following expression: 



V ~ \ co cos 



HI- 1 ) 1 -Hit!) 



cos 



(6) 



Here, we defined the polarization P = (D l F - D\)l(D l F + D F ) 
and assumed < P ^ 1 (large polarization). We see that in 
this limit the voltage increases roughly linearly with co, consis- 
tent with left panel of Fig.|2| until the frequency reaches a crit- 
ical value co cl = 2[(\+P) 2 +(\-P) 2 cos 6)] I [AP+(\-Pf sin 2 6>] 
corresponding to the condition \co sin 2 (6/2) + V\ = 1. At this 
frequency, the term 1^ starts contributing to the total current /, 
and the voltage starts decreasing monotonically with increas- 
ing co, as depicted in Fig. [2] The critical voltage in the circuit 
at co = co c \ reaches a value V c \ given by 



Vc 



4Pcos6 



cl 



4/> + (l-/>) 2 sin 2 



(7) 



which approaches unity as 6 — > 0. This is the maximum 
voltage achievable in the circuit, which, in physical units, is 
bounded by the superconducting gap As/2. Note that V gener- 
ally increases with decreasing angle 6, as shown in Fig.[2| be- 
ing in stark contrast to the typical F|I|F magnetic junction with 
one ferromagnet being free and one pinned, where V oc sin 2 6, 
i.e., vanishing with microwave power, as the precession angle 
# — > 0. Moreover, in the present setup, the voltage does not 
depend on the orientation of the precession axis, as opposed 
again to the F|I|F, where the induced voltage is sensitive to the 
relative orientation of the axis of precession with respect to 
the pinned reference ferromagnet 00. 

For co > co c i, we find that V decreases with increasing co, 
and it is given approximatively by the following expression: 



V 



1 + 



2P , 6\ , 6 

co sin - - co sin -, 

1+P 2 2 



(8) 



where we neglected corrections of the order of 
I / co cos 2 (6/2) <<c 1. Moreover, co sin 2 (6/2) ~ 1, since, in this 
regime, \co sin 2 (6/2) + V\ > 1 and \co sin 2 (6/2) - V\ < 1, while 
< V < 1 for all frequencies co. 

As the frequency increases further beyond co c i, the volt- 
age continues to decrease according to Eq. ^ until the term 
III m me tota l current / is activated, which happens when 
|6i;sin 2 (6/2) - V\ = 1. This determines the second transition 
frequency, co c2 « (1 + P) 2 /(\ + 2P)sm 2 (6/2), at which the 
voltage is given by 



V c2 



2P+V 



(9) 



neglecting corrections of order 1 /co cos 2 (6/2). According to 
Eq. ([9]), the maximum voltage at this transition point is V C 2 ~ 
1/3, corresponding to P = 1. 

Finally, for co > co C 2, all terms in Eq. §5§ contribute to the 
charge current /, and the voltage tends to zero as co — > oo. 
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FIG. 3. Dependence of the induced voltage V on oj (left panel), 
and cos 6 (right panel) for different temperatures. The increasingly 
lighter gray curves correspond to T/T c = 0.1,0.3,0.5,0.7, and 0.9 
with cos0 = 0.8 (left) and oj = 3 (right). Here c\ and cl label the 
critical points at which new charge transport channels are activated 
[see Eq. and Eq. (9)] at T = 0. In these plots, we set P = 2/3, 
and expressed V and oj in units of A s (0)/2. 
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FIG. 4. Dependence of the induced voltage V on temperature for 
different values of oj (left panel) and angle 6 (right panel). Left: 
Angle dependence of the induced voltage, for oj - 1.1, 2, 3, 5, and 10, 
lighter gray corresponding to increasing frequency, with cos 6 - 0.9. 
Right: The increasingly lighter gray curves correspond to cos 6 - 
0.9, 0.7, 0.5, 0.3, and 0.2, with oj = 3. In these plots, we set P = 2/3 
and expressed V and oj in units of A s (0)/2. 



Specifically, for oj » oj C 2, we find the following approxima- 
tion for the voltage: 



V: 



oj sin 2 (0/2) 



(10) 



The vanishingly small voltage at large oj reproduces the nil re- 
sult of ferromagnet|normal-metal junctions |8], since the the 
superconducting correlations become unimportant at frequen- 
cies oj » As . 

Note that at each transition point from one regime to an- 
other, the slope of V with respect to oj is discontinuous, as 
opposed to the case of the normal magnetic junctions, where 
it is constant for typical microwave frequencies oj <^c A/7 l8l. 
The resulting voltages are limited by the driving frequency, 
achieving values ~ oj when oj ~ 1 (or oj ~ A s /2, in phys- 
ical units) which are well within the experimental reach for 
typical superconducting gaps on the order of a few Kelvins, 



corresponding to a fraction of an meV. This is significantly 
larger than ~ 1 yueV voltages induced by magnetic dynamics 
in normal metals |7'|. 

At finite temperatures, we expect rounding off of the sharp 
features present at the transitions between different aforemen- 
tioned regions as well as a reduction of the induced voltage, 
due to the diminished superconductivity. In order to calcu- 
late explicitly the pumped current at T ^ using Eq. ([4]), we 
need to account not only for the thermally broadened Fermi 
distribution but also for the T dependence of the gap A^ (T), 
which closes at T c . We write A s (T) = A s (0)F(T), where F(T) 
is a dimensionless function, which can be found numerically 
from the self-consistent gap equation in the BCS theory fl2l . 
In Fig. [3] we plot the corresponding dependence of V on oj 
(left) and cos (right) at different temperatures, with the black 
curves showing the T = result. We see that the signal be- 
comes visibly reduced as the temperature increases from zero 
and, moreover, it singularly changes the cos — > 1 behavior 
such that V —> for any T ± instead of a finite V —> 1 at 
precisely T = 0. 

In the left panel of Fig. [4j we plot the induced voltage as 
a function of T for different frequencies oj (left) and cos0 
(right). As expected, all sharp features are smoothed out by 
finite T, with the signal eventually vanishing as T — > T c . At 
subcritical temperatures, the voltage shows a nonmonotonic 
behavior as a function of T, for a wide range of frequencies 
and angles, surprisingly reaching values in excess of the T = 
result. This is attributed to thermal activation of otherwise 
closed pumping channels. Had the gap been the same at all 
temperatures, the voltage would have exhibited an even more 
dramatic increase as the temperature is raised from T = 0, 
which is suppressed by the thermal reduction of the gap. 

In conclusion, we analyzed the dynamics of a ferromagnet 
tunnel-coupled to a conventional superconductor. We find a 
large (compared to the normal junction) nonadiabatic electri- 
cal polarization induced in an open circuit, when oj ~ A^ . We 
speculate that it could be useful, in practice, as a local probe 
for magnetic dynamics. This voltage can be easily understood 
to stem from the fact that the superconductor effectively be- 
haves as a static reference ferromagnet in the rotating frame, 
and we predict the maximum induced voltage of order of A^ , 
for small precession angles and temperatures T <T C . 
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